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Abstract: In this paper, a matrix representation for the differential resultant of two generic 
ordinary differential polynomials /i and /2 in the differential indeterminate y with order 
one and arbitrary degree is given. That is, a non-singular matrix is constructed such that 
its determinant contains the differential resultant as a factor. Furthermore, the algebraic 
sparse resultant of /i, /2, <5/i, <5/2 treated as polynomials in y, y', y" is shown to be a non-zero 
multiple of the differential resultant of /i,/2- Although very special, this seems to be the 
first matrix representation for a class of nonlinear generic differential polynomials. 
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1 Introduction 

Multivariate resultant, which gives a necessary condition for a set of n -|- 1 polynomials in 
n variables to have common solutions, is an important tool in elimination theory. One of 
the major issues in the resultant theory is to give a matrix representation for the resul- 
tant, which allows fast computation of the resultant using existing methods of determinant 
computation. By a matrix representation of the resultant, we mean a non-singular square 
matrix whose determinant contains the resultant as a factor. There exist stronger forms 
of matrix representations. For instance, in the case of two univariate polynomials in one 
variable, there exist matrix formulae named after Sylvester and Bezout, whose determinants 
equal the resultant. Unfortunately, such determinant formulae do not generally exist for 
multivariate resultants. Macaulay showed that the multivariate resultant can be represented 
as a ratio of two determinants of certain Macaulay matrixes [15]. D'Andrea established a 
similar result for the multivariate sparse resultant [7] based on the pioneering work on sparse 
resultant [2t ll0t [T6]. This paper will study matrix representations for differential resultants. 

Using the analogue between ordinary differential operators and univariate polynomi- 
als, the differential resultant for two linear ordinary differential operators was studied by 
Berkovich and Tsirulik [l] using Sylvester style matrices. The subresultant theory was first 
studied by Chardin [5j for two differential operators and then by Li [H] and Hong [11] for 
the more general Ore polynomials. 

For nonlinear differential polynomials, the differential resultant is more difficult to define 
and study. The differential resultant for two nonlinear differential polynomials in one variable 
was defined by Ritt in [181 p. 47]. In [22\ p. 46], Zwillinger proposed to define the differential 
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resultant of two differential polynomials as the determinant of a matrix following the idea 
of algebraic multivariate resultant, but did not give details. General differential resultants 
were defined by Carra Ferro using Macaulay's definition of algebraic resultants [1]. But, the 
treatment in [Jj is not complete, as will be shown in Section 2.2 of this paper. In [21], Yang, 
Zeng, and Zhang used the idea of algebraic Dixon resultant to compute the differential 
resultant. Although very efficient, this approach is not complete and does not provide a 
matrix representation for the differential resultant. Differential resultants for linear ordinary 
differential polynomials were studied by Rueda-Sendra [20]. In [19], Rueda gave a matrix 
representation for a generic sparse linear system. In [8], the first rigorous definition for 
the differential resultant of n + 1 differential polynomials in n differential indeterminates 
was given and its properties were proved. In [121113] . the sparse resultant for differential 
Laurent polynomials was defined and a single exponential time algorithm to compute the 
sparse resultant was given. Note that an ideal approach is used in [8lll21[T3]. and whether 
the multivariate differential resultant admits a matrix representation is left as an open issue. 

In this paper, based on the idea of algebraic sparse resultants and Macaulay resultants, 
a matrix representation for the differential resultant of two generic ordinary differential 
polynomials /i , /2 in the differential indeterminate y with order one and arbitrary degree 
is given. The constructed square matrix has entries equal to the coefficients of /i,/2, their 
derivatives, or zero, whose determinant is a nonzero multiple of the differential resultant. 
Furthermore, we prove that the sparse resultant of fi, f2,^fi,^f2 treated as polynomials in 
y,y',y" is not zero and contains the differential resultant of /i,/2 as a factor. Although 
very special, this seems to be the first matrix representation for a class of nonlinear generic 
differential polynomials. 

The rest of the paper is organized as follows. In Section 2, the method of Carra Ferro 
is briefly introduced and the differential resultant is defined following [8]. In Section 3, a 
matrix representation for the differential resultant of two differential polynomials with order 
one and arbitrary degree is given. In Section 4, it is shown that the differential resultant can 
be computed as a factor of a special algebraic sparse resultant. In Section 5, the conclusion 
is presented and a conjecture is proposed. 

2 Preliminaries 

To motivate what we do, we first briefly recall Carra Ferro's definition for differential resultant 
and then give a counter example to show the incompleteness of Carra Ferro's method when 
dealing with nonlinear generic differential polynomials. Finally, definition for differential 
resultant given in [8] is introduced. 

2.1 A sketch of Carra Ferro's definition 

Let K be an ordinary differential field of characteristic zero with 5 as a derivation operator. 
K{y} = K[5'^y, n € N] is the differential ring of the differential polynomials in the differential 
indeterminate y with coefficients in K. Let pi (respectively P2) be a differential polynomial 
of order m and degree di (respectively of order n and degree d2 ) in K{y}. According to 
Carra Ferro [1], the differential resultant oi pi,p2, denoted by 6R{pi,p2), is defined to be the 
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Macaulay's algebraic resultant of m + n + 2 differential polynomials 

V{pi ,P2) = {5yi,6^~^pi , . . . , Pi , O2, S"'~'p2, ...,P2} 

in the polynomial ring Sm+n = K[y, 6y, . . . , 6'"^~^""y] in m + n + 1 variables. 
Specifically, let 



Z) = l + (n+l)((ii-l) + (m+l)((i2-l), L 



m + n + I + D 
m + n + 1 



Let Ui = 6''y for all i = 0,1,... ,m + n. For each a = (ao, . . . , a^+n) € N"*+"+\ Y"' = 
Vo^ ' ' ' Vm+n is a power product in Sm+n- M^j^^j^^ stands for the set of all power products 
in Sm+n of degree less than or equal to D. Obviously, the cardinality of M^_,_„_,_]^ equals L. In 
a similar way it is possible to define M^^^^^^ which has Li = i^'^^^^^^'^^) monomials, and 
M^^^^i which has L2 = ("^^+^^_,^"'^') monomials. The monomials in M^^^.^^, M^^^^^^^ 
and Mj^^^^]^ are totally ordered using first the degree and then the lexicographic order 
derived from yo < yi < • • • < Vm+n- 

Definition 2.1 The {{{n + l)Li + (m + 1)^2) x L)-matrix 

M{6,n,m)= M(5"pi ,...,5pi,pi, 6'^p2 , . . . ,6p2,P2), 

is defined in the following way: for each i such that (j — l)Li < i < jLi the coefficients of 
the polynomial y"(5""^ ■'pi are the entries of the i-th row for each V £ M^~^]^^ and each 
j = 1, . . . , n + 1, while for each i such that 



(n + l)Li + (j - n - 2)L2 < i < (n + l)Li + {j - n - 1)L 
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the coefficients of the polynomial Y°'5^^^^'^~^p2 are the entries of the i-th row for each 
Y"" G M^^^^^ and each j = n + 2, . . . ,m -\- n + 2, that are written with respect to the power 
products in M^_^_^_^_-^ in decreasing order. 

Definition 2.2 The differential resultant of pi and p2 is defined to be 

gcd(det(P) : P is an (L X L)-suhmatrix of M{5,n,m)). 



2.2 A counter example and definition of diff"erential resultant 

In this subsection, we use Carra Ferro's method to construct matrix formula of two 
nonlinear generic ordinary differential polynomials with order one and degree two, 

91 = aoyf + aiyiy + 02^^ + asyi + a^y + 05, 

92 = boyf + biyiy + 62y^ + ^'3^1 + hv + ^5, (1) 

where, hereinafter, yi = 6y and Oj, bi with i = 0, . . . ,5 are generic differential indeterminates. 

For differential polynomials in ([1]), we have di = d2 = 2,m = n = 1, D = 5, L = 56, and 
Li = L2 = 20. The set of column monomials is 

Ml = {yl,ylBlylBlylBly2Bl,Bl}, 
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where, and throughout the paper, B = yi,y2} and Bf denotes all monomials of total 
degree less than or equal to j in the first i elements of B. For example, B2 = {l,y,y'^} and 
i?| = {l,y,yi,y'^,yyi,yl}. Note that the monomials of are M3 = {yf, 2/2-^3 > 2/2 -Bf} = 



R3 



According to Definition 12. H M{6, 1, 1) is an 80 x 56 matrix 
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Obviously, the entries of the first column are all zero in M{5, 1, 1), since the monomial 
never appears in any row polynomial Y * f, where the monomial Y € M3 and / € 

{5gi,gi, (5(72) 52}- Consequently, the differential resultant of gi,g2 is identically zero according 

to this definition. 

Actually, the differential resultant is defined using an ideal approach for two generic 
differential polynomials in one differential indeterminate in [18] and re + 1 generic differential 
polynomials in n differential indeterminates in [8]. / is said to be a generic differential 
polynomial in differential indeterminates Y = {yi, . . . ,y„,} with order s and degree m ii f 
contains all the monomials of degree up to rei, in yi , . . . , ?/„ and their derivatives up to order 
■s. Furthermore, the coefficients of / are also differential indeterminates. For instance, gi 
and g2 in ([T]) are two generic differential polynomials. 

Theorem 2.3 ( [8]) Let pQ,pi, . . . ,pn be generic differential polynomials with order Si and 
coefficient sets Uj respectively. Then [po,Pi, ■ ■ ■ ,Pn] is a prime differential ideal in Q{Y, uq, 
. . . , u„}. And 

[po,Pi, ■ ■ ■ ,Pn] n Q{uo, . . . ,u„} = sat(R(uo, . . . ,u„)). (2) 

is a prime differential ideal of codimension one, where R is defined to be the differential 
sparse resultant of po,pi, . . . ,pn, which has the following properties 

a) R(uo, ui, . . . , u„) is an irreducible polynomial and differentially homogeneous in each 

Uj. 

b) R(uo, ui, . . . , u„) is of order hi = s — si in Wi (z = 0, . . . , n) with s = X^"=o ^i- 
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c) The differential resultant can be written as a linear combination of pi and the deriva- 
tives of Pi up to the order s — Si. Precisely, we have 

n s—Si 

R(uo, ui, . . . , u„) = ^ ^ hijp[^^ 

i=0 j=Q 

where hij € Q[Y, Y^^) . . . , Y(^), u[,'"'°\ . . . , ui'~'"^]. 

3 Matrix formula for differential polynomials 

In this section, we will give a matrix representation for the following generic differential 
polynomials in y 

fi = aydivt + ayfi-\yyt~^y H ^ ao, (3) 

/2 = byd^yf'' + byd^-i^yf^'^y H h 60, 

where yi = 5y, I < di < d2, and a , . . . ,ao,b da , • • • , 60 are differential indeterminates. 
3.1 Matrix construction 

In this subsection, we show that when choosing a proper column monomial set, a square 
matrix can be constructed following Macaulay's idea |15| . 

By c) of Theorem 12.31 the differential resultant for /i,/2 can be written as a linear 
combination of /i , /2 , 5/i , 5/2 which are treated as polynomials in variables y, 2/1,2/2 = ^yi- 
So, we will try to construct a matrix representation for the differential resultant of /i,/2 
from these four polynomials. 

From Section 2, it is easy to see that the main problem for Carra Ferro's definition is 
that the matrix M{5, n, m) contains too many columns. Or equivalently, the monomial set 
-^m+n+i used to represent the columns is too large. 

Consider the monomial set 

E = B^Uy2Bt^ (4) 

with D = 2di + 2d2 — 3. We will show that if using S as the column monomial set, a 
nonsingular square matrix can be constructed. 

Define the main monomial of polynomials pi = 6fi,p2 = 5/2, P3 = /i,P4 = /2 to be 

mm(pi) = y2?/^~\nim(p2) = yf'^,m.m{ps) = y'^\mm(p4) = 1. (5) 

Then, we can divide £ into four mutually disjoint sets £" = 5i U ^2 U ^3 U ^4, where 

5i = {y" G f : mm(pi) divides F"}, 

^2 = G £ : mm(pi) does not divide but mm(p2) does}, (6) 

53 = {y" € £ : mm(pi), mm(p2) do not divide y" but mm(p3) does}, 

54 = {y" G £ : mm(pi), mm(p2), mm(P3) do not divide y"}. 
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As a consequence, we can write down a system of equations: 



y°/mm(pi) *pi = 0, for G 5i, 

y°/mm(p2)*P2 = 0, for G ^2, 

y°/mm(p3)*P3 = 0, for G ^3, (7) 

y"/mm(p4) * P4 = 0, for G 54. 

Observe that the total number of equations is the number of elements in S and denoted by 
N = (D + l)'^. 

Regarding the monomials in ([7]) as unknowns, we obtain a system of N linear equations 
about these monomial unknowns. Denote the coefficient matrix of the system of linear 
equations dTj) by -Ddi,d2 whose elements are zero or the coefficients of fi and 6fi,i = 1,2. 

Note that the main monomials of the polynomials are not the maximal monomials in the 
sense of Macaulay [15], so the monomials on the left hand side of (jT]) may not be contained 
in £. Next, we prove that this does not occur for our main monomials. 

Lemma 3.1 The coefficient matrix D(i^^d2 of system ^ is square. 

Proof: The coefficient matrix of system ([7]) has N = \S\ rows. In order to prove the lemma, it 
suffices to demonstrate that, for each G 5^, i = 1, . . . , 4, all monomials in [y"/mm(pj)] *pi 
are contained in £. Recall that £ = U y2B^~^. Then by ([6]), one has 

51 = B^-"^^ * mm(pi) = 5f "'^1 * mm((^/i), 

52 = ^ mm(p2) = B^-^'' * mm{6f2), 

53 = Ti * mm(p3) = Ti * mm(/i), (8) 

54 = T2 * mm(p4) =T2* mm(/2). 



where 



d2-l 

jD-di-l-i 



^1 = ( U ylBt''-') U [y^ U yl^i 

I i=0 i=0 

T2 = \([}y{Bt~') U U^1^2^") \ ■ (9) 

I i=0 i=0 J 

Note that the representation of ^2 in ^ is obtained with the help of the condition di < ^2- 
Hence, the equations ([7]) become 

B^-^^ * 5fi = 0, 

B^-"'' * 5f2 = 0, (10) 
Ti * /i = 0, 
T2 * /2 = 0. 

Since the monomial set of 5fi is -63^ U 1/2 * -^3^"^) t^is monomial set of B^~'^^ * 5fi is 
B^~'^^ * U y2 * -63^"^) = -Bf U y2-Bf = £. So monomials in the first set of equations 
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in (llOp are in 8. Since the monomial set of /i is .63^ = Uj^Qy^-B2^ > the monomial set of 
Ti * h is Tn U 2/2ri2, where Tn = Uk=o'~^ ViB?'"' and Tia = U'lV' yi^2''"''"'- Since 
di > 1 and d2 > 1, we have di + d2 — ^ < D = 2di + 2d2 — 3 and hence Tn C . Since 
di > 1 and d2 > 1, we have 2(ii - 2 < D - 1 = 2di + 2^2 - 4 and hence T12 C B^^^. As a 
consequence, Tn U 2/2 7i2 C f and the monomials in the third set of equations in (jlOp are in 
i5. Other cases can be proved similarly. Thus all monomials in the left hand side of (jlOp are 
in £. This proves the lemma. □ 

It is worthy to say that, due to the decrease of the number of monomials in £ compared 
with the method by Carra Ferro, the size of the matrix Dci^^(i2 decreases significantly. 

3.2 Matrix representation for differential resultant 

In this section, we show that det(D^j^rf2) is not identically equal to zero and contains the 
differential resultant factor. 

Lemma 3.2 det{D^^ (i2) identically equal to zero. 

Proof: It suffices to show that there exists a unique monomial in the sense that it is 
different from all other monomials in the expansion of (iel{D(i^^d2)- 

The coefficients of the main monomials in 5fi , 5f2 , /i , /2 are respectively 

the coefficient of mm(5/i) = y2yi^~^, 

the coefficient of mm((5/2) = yf'^, (11) 

the coefficient of mm(/i) = y'^^ , 
the coefficient of mm(/2) = 1. 

We will show that the monomial (a di)^^{6b ti2)"'^(%rfi)"^(^o)""' is a unique one by the 
following four steps, where rij is the number of elements in Si with i = 1, . . . , 4. From (jlOp . 
m = 1^3^-'^^, n2 = \B^-'''\, ng = iTil, n4 = \T2\. 

1. Observe that, in dfi, Uyd^ only occurs in the coefficient of yiy'^^"^ with the form 
Uyd-^ + SOy^yd^^i. Fu r t h 6 r m o r e , SOy^yd^-i only occurs in this term given by 6fi and no other 
places of Dd^ ci2- So using the transformation 

^^yiy'^i-^ ~ Sij/'*!"^ ~ diQyd^, with other coefficients unchanged, (12) 

where Cy^yd^-i is a new differential indeterminate, Dd^^d2 is transformed to a new matrix 
which is singular if and only if the original one is singular. Still denote the matrix by Z)^^ ,^2- 
From (jlOp . for a monomial M G Ti, ayd^ is the coefficient of the monomial My'^^ in 
each polynomial Ti * fi and hence in each corresponding row of 1)^1,^2 • Then Uydj is in 
different rows and columns of D^^^^j, and this gives the factor (Oj^di)"^. Delete those rows 

and columns of D^^^d^ containing Qyd^ and denote the remaining matrix by D^^^^. From 

(jlOp . the columns deleted are represented by monomials y'^^Ti. So, D^^^^^ i^ ®*i^^ ^ square 
matrix. 





a d, 
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5b dn +h dn 
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bo 
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2. Let M G bP^'^K The term a di occurs in M * 6fi as the coefficient of the monomial 



Vl 

y2yi^~^ M , or equivalently it occurs in the columns represented by ?/2?/f^~^M. This gives the 
factor (a di . It is easy to check that a does not occur in other places of d'^_^^^ . From the 

definition for Ti Q, the columns deleted in case 1 correspond those columns represented by 
monomials of the form y2^yi^y^^ where either k2 = and ki < d2 or k2 = 1 and ki < di — 1. 
Then {y'^^Ti} n {y2yf^~^ -<^i| — 0^ qj. equivalently those columns of D^-^ ^^ containing 

a di are still in D^}\ . Similar to case 1, one can delete those rows and columns of D^P, 

y^J- ai,a2 ' fil>"2 

(2) 

containing a^d-^ and denote the remaining matrix by D\^^ which is still a square matrix. 
From (jlOp . the columns deleted are represented by monomials y2yi^~^ B^~'^^ . 

(2) 

3. At the moment, -D^^^^j o^^^Y contains coefficients of /2 and 6/2- Observe that bo only 
occurs in the rows corresponding to T2* f2, where T2 is defined in ([9]). Note that 5bo instead 
of 60 occurs in 6/2- Since {y'^^Ti} n r2 = and {y2yi^~^B^~'^^} n r2 = 0, the columns 
of -D(ii,d2 containing bQ, represented by T2, are not deleted in case 1 and case 2. Then, we 

f2) 

have the factor (^g)""'- Similarly, delete those rows and columns of -D^^^^ containing 69 and 

denote the remaining matrix by D^_^^^ which is still a square matrix. From ()10p , the columns 
deleted are represented by monomials T2. 

4. From (jlOp . the rows of D^f^^^ are from coefficients of B^~'^'^ *5/2- The term 5b d2 is in 

the coefficient of the monomial M * yf^ in M * 6f2 for M € bP~'^''^ , and 6b d2 does not occur 

in other places of M *6f2- Furthermore, since {y'^^Ti}r\{yf'^ b!^~'^^} = 0, {y2yi^~^ B^~'^^} n 
{yi^ bP'"^^} = 0, and T2 n {Vi^ bP~'^^}} = 0, the columns containing the term 6b da are not 

Vl 

deleted in the first three cases. Then, we have the factor {6b^d2 • 

Following the above procedures step by step, the coefficients of choosing main monomials 
of the polynomials /i, /2, 6fi, 6/2 occur in each row and each column of i5(ii,d2 ^iid only once, 
and the monomial {o-ydi)^^ {6b^d2)^^ {ttydi)"'^ [bo)^'^ is a unique one in the expansion of the 
determinant of D^^ ^i^. So the lemma follows. □ 

Note that the selection of main monomials in above algorithm is not unique, thus there 
may exist other ways to construct matrix formula for system ([3]). 

Corollary 3.3 Following the above notations, for any G Si, if all monomials of[Y°'/mm{pj 
Pj are contained in £{j 7^ i), then the rearranged matrix, which is obtained by replacing the 
row polynomials [Y°'/mm{pi)] * pi by [y"/mm(pj)] *Pj, is not identically equal to zero. 

Corollary 13.31 follows from the fact that the proof of Lemma 13.21 is independent of the 
number of elements in Si as long as the main monomials are the same. 

The relation between det(-Dd^.d2) and differential resultant of /i,/2, denoted by R, is 
stated as the following theorem. 

Theorem 3.4 det^Ddj^^^^) is a nonzero multiple o/R. 

Proof. From Lemma 13.21 det (1)^1,^2) nonzero. In the matrix D^-^ ,^^, multiply a column 
monomial M 7^ 1 in to the corresponding column and add the result to the constant 
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column corresponding to the monomial 1. Then the constant column becomes Y"' * pi with 
Pi = <^/i)P2 = Sf2,P3 = /i,P4 = /2 and £ Si/mm{pi), i = 1, . . . ,4. Since a determinant 
is multi-linear on the columns, expanding the matrix by the constant column, we obtain 

det{Dd,4,) = hifi + h2 6fi + /i3 /2 + /i4 Sf2, (13) 

where hj are differential polynomials. From ([2]), det{Dii^ [i^) G sat(R). On the other hand, 
from Theorem 12.31 R is irreducible and the order of R about the coefficients of /i, /2 is one. 
Therefore, R must divide det{Ddj^d2)- ^ 

From Theorem 13. 4^ we can easily deduce a degree bound = 4((ii + d2 — 1)^ for the 
differential resultant of /i and f2- The main advantage to represent the differential resultant 
as a factor of the determinant of a matrix is that we can use fast algorithms of matrix 
computation to compute the differential resultant as did in the algebraic case [3]. 

Suppose that det{Dd^^d2) is expanded as a polynomial. Then the differential resultant 
can be found by the following result. 

Corollary 3.5 Suppose det(Drfj^d2) = Y\i=i irreducible factorization of det[Dd^^d2) 

in Q[C/j,C/2], where Cf^,i = 1,2 are the sets of coefficients of fi. Then there exists a unique 
factor, say Pi, which is in [/i,/2] and is the differential resultant of fi and f2- 

Proof. From c) of Theorem 12.31 and Theorem 13.41 R € [/i,/2] and is an irreducible factor 
of dei{D(i^^(i^). Suppose dei{D(i^^(i^) contains another factor, say P2, which is also in [/i,/2]. 
Then P2 € sat(R) by ([2]). Since R is irreducible with order one and P2 is of order no more 
than one, P2 must equal R, which contradicts to the hypothesis. □ 

3.3 Example ([1]) revisited 

In this section, we apply the method just proposed to construct a matrix representation for 
the differential resultant of the system ([1]) . 

Following the method given in the proceeding section, for system ([T]), we have D = 2di + 
2d2 — 3 = 5 and select the main monomials of 6gi,6g2,gi,g2 are y2yi,yi,y'^ , 1 respectively. 
Then £ = 7/2^3 U is divided into the following four disjoint sets 

51 = y2yiBi, 

52 = yfBl 

53 = y^ [Bl U yiBl U 2/2^2'] , (14) 

54 = BluyiBluy2Bl 

Using ([7]) and regarding the monomials in £ as variables, we obtain the matrix 1)2,2 ; 
which is a 36 X 36 square matrix in the following form. 
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Bl5g2< 
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As shown in the proof of Lemma 13.21 {aoY^ {a2Y^ {b^)^ {^h^Y^ is a unique monomial in 
the expansion of the determinant of -D2,2- Hence, the differential resultant of gi and 52 is a 
factor of det(L'2,2)- Note that in Carra Ferro's construction for gi,g2, M{5, 1, 1) is an 80 x 56 
matrix, which is larger than 1)2,2 • 

In particular, suppose uq = bo = 1 and aj,6j are differential constants, i.e, 6ai = 5bi = 0, 
i = 1, ... ,5. Then 1)2,2 can be expanded as a polynomial and the differential resultant of 
91,92 can be found with Corollary 13.51 which is a polynomial of degree 12 and contains 3210 
terms. This is the same as the result obtained in 1211. 



4 Differential resultant as the algebraic sparse resultant 

In this section, we show that differential resultant of fi and /2 is a factor of the algebraic 
sparse resultant of the system {/i, /2, (^/i, '5/2}• 
4.1 Results about algebraic sparse resultant 

In this subsection, notions of algebraic sparse resultants are introduced. Detailed can be 
found in [SlElIIOlIIE] . 

A set S in R"- is said to be convex if it contains the line segment connecting any two 
points in S. If a set is not itself convex, its convex hull is the smallest convex set containing 
it and denoted by Conv(5). A set F = {ai, . . . , Om} is called a vertex set of a convex set Q 
if each point q & Q can be expressed as 

m m 

q = AjOj, with Xj = 1 and Xj > 0, 
i=i i=i 

and each aj is called a vertex of Q. 

Consider n + 1 generic sps^rse polynomials in the ctlgebraic indeterminfites x^, . . . ^Xn- 

Pi = Uio + UiiMil H h Uii.Mii., i = 1, . . . ,n + 1, 
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where Uij are indeterminates and Mj^ = ns=i ' monomials in Q[xi,...,x„] with 
exponent vectors aik = (e*'^^ , . . . , e*'^") G Z". Note that we assume each pi contains a 
constant term u^q- For a = (e^j-.-je") G Z", the corresponding monomial is denoted as 

M(a) = n:=i<. 

The finite set Ai C of all monomial exponents appearing in pi is called the support of 
Pi, denoted by supp(pi). Its cardinality is li = \Ai\. The Newton polytope Qi C R" of pi is 
the convex hull of ^j, denoted by Qj = Conv(^j). Since Qi is the convex hull for a finite set 
of points, it must have a vertex set. For simplicity, we assume that each Ai is of dimension 
n as did in [101 P-252]. Let u be the set of coefficients of pj, i = 0, . . . , n. Then, the ideal 

(pi,P2, . . . ,Pn+i) n Q[u] = (7^(u)) (15) 

is principal and the generator TZ is defined to be the sparse resultant of pi, . . . ,Pn+i jlOj 
p. 252]. When the coefficients u of pi are specialized to certain values v, the sparse resultant 
for the specialized polynomials is defined to be ^(v). The matrix representation of TZ is 
associated with the decomposition of the Minkowski sum of the Newton polytopes Qi. 
The Minkowski sum of the convex polytopes Qi 

Q = Qi-\ h Qn+i = {qi-\ h qn+i\qi e Qi}- 

is still convex and of dimension n. 

Choose sufficiently small numbers 5i > and let 5 = {6i, . . . , (5„) G R" be a perturbed 
vector. Then the points which lie in the interior of the perturbed district £" = n (Q + 5) 
are chosen as the column monomial set [3] to construct the matrix for the sparse resultant. 

Choose n + 1 sufficiently generic linear lifting functions h, ■ ■ ■ , In+i £ Z[xi, . . . , x.„] and 
define the lifted Newton polytopes Qi = = {qi,li{qi)) ■ qi G Qi} C M""^^. Let 

n+1 

Q = Y^Q,C 

i=l 

which is an (n + l)-dimensional convex polytope. The lower envelope of Q with respect 
to vector (0, . . . , 0, 1) G M"+^ is the union of all the n-dimensional faces of Q, whose inner 
normal vector has positive last component. 

Let TT : (gi, . . . ,qn+i) ^ {qi, ■ ■ ■ ,qn) be the projection to the first n coordinates from 
M""*"^ to M". Then vr is a one to one map between the lower envelope of Q and Q [3]. The 
genericity requirements on li assure that every point qon the lower envelope can be uniquely 
expressed as q = qi + ■ ■ ■ + §n+i with qi (z Qi, such that the sum of the projections under vr 
of these points leads to a unique sum of q = qi + • • • + qn+i G Q C with qi £ Qi, which 
is called the optimal (Minkowski) sum of q. For Ti C Qi, R = Y17=i -^i called an optimal 
sum, if each element of it! can be written as a unique optimal sum X]"=QPi for pi G J-^. 

A polyhedral subdivision of an ?i-dimensional polytope Q consists of finitely many n- 
dimensional polytopes Ri, . . . , Rg, called the cells of the subdivision, such that Q = RiU 
■ ■ - U Rg and for i ^ j and Ri D Rj is either empty or a face of both Ri and Rj. A polyhedral 
subdivision is called a mixed subdivision if each cell Ri can be written as an optimal sum 
Ri = Y17=i^-^i^ where each Ti is a face of Qi and n = X^^i^ dim(J-'j). Furthermore, if 
Rj = J27=i -^i another cell in the subdivision, then Ri n Rj = Yll=i{'^i ^ -^i)- ^ ^^^^ 
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Ri = X^"^;^ J-i is called mixed if dim(J-i) < 1 for all i; otherwise, it is a non-mixed cell. As a 
result of n = "^^^i dim(J'j), a mixed cell has one unique vertex, which satisfies dim(J^jp) = 0, 
while a non-mixed cell has at least two vertices. 

Recall 6 = {5i,52, ■ ■ ■ , 6n), where < di < 1. If Q = iZi U • • • U iZs is a subdivision of Q, 
then 6 -\-Q = {6 + Ri)U---U {6 + Rs) is a subdivision oi 6 + Q. 

Let g € n (Q + (5) lie in the interior of a cell 6 + J-i + ---+ J-n+i of a mixed subdivision 
for Q + 6, where J-i is a face of Qi. The row content function of q is defined as the largest 
integer such that J-i is a vertex. In fact, all the vertices in the optimal sum of p can be 
selected as the row content functions. Hence, we define generalized row content functions 
(GRC for brief) as one of the integers, not necessary largest, such that J-i is a vertex. 

Suppose that we have a mixed subdivision of Q. With a fixed GRC, a sparse resultant 
matrix can be constructed as follows. For each i = l,...,?7- + l, define the subset Si of £ as 
follows: 

Si = {qe£\GRC{q) = {i,jo)}, (16) 
where jo € {1, . . . ,mj}, mj is the number of vertexes of Qi, and we obtain a disjoint union 
for £: 



£ = SiU---USn+i. (17) 

For q £ Si, let q = qi + ■■■ + q-n+i G Q be an optimal sum of q. Then, q^ is a vertex of Qi 
and the corresponding monomial M{qi) is called the main monomial of pi and denoted by 
mm(pj), similar to what we did in Section 3. Main monomials have the following important 
property [6l p350]. 

Lemma 4.1 If q £ Si, then the monomials in {M{q)/mm{pi))pi are contained in £. 

Now consider the following equation systems 

{M{q)/mm{pi))pi, q e Si,i = 1, . . . ,n + 1. (18) 

Treating the monomials in £ as variables, by Lemma 14. H the coefficient matrix for the 
equations in (jlSp is an \£\ x \£\ square matrix, called the sparse resultant matrix. The sparse 
resultant oi pi,i = 1, ... ,n + 1 is a factor of the determinant of this matrix. 

In [2l[3], Canny and Emiris used linear programming algorithms to find the row content 
functions and to construct Si. We briefly describe this procedure below. 

Now assume Qi has the vertex set Vi = {an, . . . , aj m, }• A point q £ ri{Q + 6) implies 
that q £ a + 6 with a cell a £ Q. In order to obtain the generalized row content functions 
of q, we wish to find the optimal sum of g' — 5 in terms of the vertexes in Vi. Introducing 
variables Xij,i = 1, . . . ,n + 1, j = 1, . . . ,mi, one has 

n+1 n+1 rrii mi 

q — 6 = qi = Xij Oij, with Ajj = 1 and Ajj > 0. (19) 

i=l 1=1 j=l j=l 

On the other hand, in order to make the lifted points lie on the lower envelope of Q, one 
must force the "height" of the listed points minimal, thus requiring to find Xij such that 

n+1 mi 

Ajj (oj j ) to be minimized (20) 

i=i j=i 
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under the linear constraint conditions ()19p . where li{aij) is a random Hnear function in aij. 
For q £, let A^^- be an optimal solution for the linear programming problem (j20p . Then 

q — 6 = J27=i Qi where q* = X*jli{aij). a*j^ is a vertex of Qi if and only if there exists 

a jo such that A^^^ = 1 and A*^ = for j ^ Jq. In this case, the generalized row content 
function of q is (i,jo) and mm(pj) is M{a*j^). It is shown that when the lift functions /j are 
general enough, all Si can be computed in the above way [3]. 

In order to study the linear programming problem ()20p . we need to recall a lemma about 
the optimality criterion for the general linear programming problem 

T 

mm z = c X 

X 

subject to Ax = b, with I < x < u, (21) 

where ^ is an m x n rectangular matrix, 6 is a column vector of dimension m, c and x are 
column vectors of dimension n, and the superscript T stands for transpose. In order for 
the linear programming problem to be meaningful, the row rank of A must be less than the 
column rank of A. We thus can assume A to be row full rank. Let 77,1, .. . ,nm be linear 
independent columns of A. Then the corresponding x^ , ■ ■ ■ , Xn^ are called basic variables 
of X. Let B be the matrix consisting of the ni, . . . , nm columns of A. Then S is an m x m 
invertible matrix. Lemma [4. 21 below gives an optimality criterion for the linear programming 
problem (f2T]) . 

Lemma 4.2 ( J^) Let xb be a basic variables set of x, where B is the corresponding co- 
efficient matrix of xb- If the corresponding basic feasible solution xb = B~^b > and the 
conditions cbB~^ A—c < hold, where cb is the row vector obtained by listing the coefficients 
of Xb in the object function, then an optimal solution for the linear programming problem 
h21\) can be given as xb = B^^b and all other Xi equals zero, which is called the optimal 
solution determined by the basic variables xb- 



4.2 Algebraic sparse resultant matrix 

In this subsection, we show that the sparse resultant for /i, /2, 5fi , 5/2 is nonzero and contains 
the differential resultant of /i and /2 as a factor. 

For the differential polynomials /i and /2 given in ([3]), consider the pi = Sfi,p2 = 
^f2,P3 = fi,P4 = /2 as algebraic polynomials in y, 2/1,2/2- The monomial sets oiSfi,Sf2, /i, /2 
are U 2/2 * ^3'"\ U 2/2 * B^"^"^, B^\ and B^'^ respectively. For convenience, we will 
not distinguish a monomial M and its exponential vector when there exists no confusion. 
Then the Newton polytopes for Sfi,6f2, fi, /2 are respectively, 

Qi = Conv(sup(5/i)) = Conv(5^i U 2/2 * ^3'"^) C Il^ 

Q2 = Conv(sup((5/2)) = Conv(B3^ U 2/2 * ^3'^^) C Il^ 

Q3 = Conv(sup(/i)) = Conv(B^i) C Il^ (22) 

Q4 = Conv(sup(/2)) = Conv(B^2) C Il^ 

The Newton polytopes Qi and are shown in Figure 1 (for di = 5) while Q2 and (54 have 
similar polytopes as Qi and but with different sizes respectively. 
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Figure 1: The Newton polytopes Qi and Q3 



Let the Minkowski sum Q = Qi + Q2 + Q3 + Qi- In order to compute the column 
monomial set, we choose a perturbed vector 5 = {61, 82, (Js) with < 5,, < 1 with i = 1, 2, 3. 
Then the points in H {Q + 6) is easily shown to be yyiy2£ where £ is given in (j4]). Note 
that using £ or yyiy2£ as the column monomial set will lead to the same matrix. 

The vertex sets of Qi, denoted by Vj, are respectively 

Vi :={(0,0,0),(0,0,l),(0,di -1,1), (0,di,0),(di -1,0,1), (di, 0,0)}, 
V2 := {(0, 0, 0), (0, 0, 1), (0, d2 -1,1), (0, d2, 0), {d2 - 1,0, 1), (da, 0, 0)}, 
V3 :={(0,0,0),(0,di,0),(di,0,0)}, 
V4 :={(0,0,0),(0,d2,0),(d2,0,0)}. 

Let the lifting functions be li = {Ln, Li2, Li^), i = 1, ... ,4, where Lij are parameters to be 
determined later. From (120p . the object function of the linear programming problem to be 
solved is 

min(Ai2Li3 + Ai3[ii2(di - 1) + L13] + \uL12d1 + Ai5[Lii(rfi - 1) + L13] + Aigiiidi 

+ X22L23 + A23[i22(d2 ~ 1) + -^23] + A24.^22'^2 + A25[i2l((i2 ^ 1) + -^23] + ^26^21^2 

+ A32-Zj32(il + A33L3iC?i (23) 

-\-Xi2Li2d2 + A43L4iC?2) 

under the constraints 

Ai = Ai5(di - 1) + Aiedi + A25(d2 - 1) + A26d2 + Aaadi + A43d2, 

A2 = Ai3(di - 1) + Ai4di + \22.{d2 - 1) + A24d2 + A32dl + A42d2, 

^3 = Ai2 + Ai3 + Ai5 + A22 + A23 + A25, (24) 

nii 

^Aij = 1, z = 1,...,4, 

J=l 

Ajj > 0, i = 1, . . . , 4, j = 1, . . . , mj with mi = m2 = 6, m^, = = 3, 
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where Ai = ei — di, A2 = 62 — 62, ^3 = £3 — S3 with (ei, £2, £3) G Z'^ fl (Q + 6). According to 
the procedure given in Section 14.11 we need to solve the hnear programming problem (j23p 
for each {£1,82,83) G D {Q + 5). Note that Lij are parameters. What we need to do is 
to show that there exist Lij such that the solutions of (j23p make the corresponding main 
monomials to be the ones selected by us in ([5]). More precisely, we need to determine Lij 
such that for each q G Z"n (Q + J), the optimal solution for the linear programming problem 
(|23p consists one of the following cases: 

A13 = 1 implies GRC((7) = (1,3), the vertex is {0,di — 1, 1), and mm((5/i) = y2yi^~^, 

A24 = 1 implies GRC((7) = (2,4), the vertex is (0,^2,0), and mm((5/2) = yf^, 

A33 = 1 implies GRC((7) = (3,3), the vertex is ((ii,0,0), and mm(/i) = y"^^, 

A41 = 1 implies GRC((7) = (4, 1), the vertex is (0,0,0), and mm(/2) = 1. 

The following lemma proves that the above statement is valid. 

Lemma 4.3 There exist Lij such that the optimal solution of the corresponding linear pro- 
gramming problem \2'J\) can he chosen such that the corresponding main monomials for 
/i, 72,5/1,5/2 are mm(/i) = y^^S mm(/2) = 1, mm{5fi) = 2/22/i'~S mm{6f2) = yf^ re- 
spectively and £ can be written as a disjoint union £ = SiU S2U S3L) S4, where Si is defined 
in ([7^. 

Proof. We write the linear programming problem as the standard form ()2ip . It is easy to 
see 

c = (0, Li3, {di - l)Li2 + -Lis, diLi2, {di - l)Lii + L13, diLn, 

0, -^^23, {d2 — 1)^^22 + -^^23, C^2-^^22, (^2 " 1)-^^21 + -^^23, ^2-^^21, 

0, diL32, diLsi, 0, ^2-^^42, (^2-^41)- 

Let 6 = {61,62,63) be a sufficiently small vector in sufficiently generic position, the validity 
of 6 is analyzed in [3]. Then 
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--di 


-1 


d2 


= d2 


-1, 


which 


is a 


7 X 18 


matrix and b 




(^1, 


A2, 


^3, 


1,1,1,1 



easy to see that the rank of A is 7, since di > 1. 

From gD, we have f = Z'^n{Q + 6) = yyiy2{B^ Uy2B^~^), where D = 2(^1 + 2^2-3. We 
will construct a disjoint union £ = 5i U ^2 U Ss U S'4 like (fT7|) such that the corresponding main 
monomials are respectively mm(5/i) = 2/22/1^""^, mm(5/2) = y^^,mm(/i) = y'^^, mm(/2) = 1. 

Four cases will be considered. 
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Case 1. We will give the conditions about Lij under which mm(J/i) = ?/22/i^~ , or 
equivalently, the linear programming problem (|23p has an optimal solution where A13 = 1. 
As a consequence, 5i will also be constructed. 

As shown by Lemma 14.21 an optimal solution for a linear programming problem can 
be uniquely determined by a set of basic variables. We will construct the required optimal 
solutions by choosing different sets of basic variables. Four sub-cases are considered. 

1.1. Selecting basic variables as vetn = {A13, A23, A24, A32, A33, A41, A43} while other 
variables are nonbasic variables and equal to zero. Due the constraint Aii+Ai2 + - • •+A16 = 1, 
we have A13 = 1. Then for any such an optimal solution of the linear programming problem 
([23]). in the optimal sum of any element q = qi + q2 + <!?, + q^, Qi = (0, di — 1, 1) is a vertex 
of Qi, mm((5/i) = 1/22/1 ^^"^ corresponding q belongs to Si as defined in (fTHl) . 

We claim that the basic feasible solutions in vetn must be nondegenerate meaning that 
all basic variables are positive, that is, xb = B~^b > 0. A mixed cell R = Ylt=i'^i^ where 
J-i is a face of Qi, must satisfy the dimension constraint "^f^i dim(J-'j) = 3. In the cell 
corresponding to the basic variables vetn, Ti = (0,(ii — 1, 1) is a vertex of Qi, T4 is a one 
dimensional face of of the form A41V41 + A43V43, where V41 = (0,0,0), V43 = ((^2,0, 0), 
and A41 + A43 = 1. 7^2 and are one dimensional faces of Q2 and Q3 respectively. In order 
for the dimension constraint X^i=i dim(J^j) = 3 to be valid, the claim must be true. For 
otherwise, one of the variables in vetn must be zero, say A41 = 0. Then A43 = 1 and 
becomes a vertex, which implies X]i=i dim(J^j) < 3, a contradiction. 

From Lemma 14. 2| the coefficient matrix of basic variables in ()24p is 
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with rank(i?ii)=7. For all (^1,^2,^3) and b = (^1,^2,^3,1,1,1,1), the requirement 
B^^b > in Lemma 14.21 gives 

1< ^3 < 2, di + (i2 < ^2 + ^3 < 2di + d2, 
2di + d2 < Ai+ A2 + A3 < 2di + 2d2. 

Substituting Ai = ei — 5i, A2 = £2 — ^2, A3 = £3 — 83 into the above inequalities and 
considering that (£1,625 £3) are integer points, we have 

£3 = 2, e2 = di + d2-l,...,2di + d2- 2, (25) 
£1 + £2 = 2(ii + (i2 - 1, . . . , 2di + 2d2 - 2. 

On the other hand, cb^^ = {{di - l)Li2 + -^^13, (c?2 - l)-^^22 + -i^23, c?2-i^22, ^1-^^32, di-Z>3i, 0, 
<^2-^4i)- After simplification and rearrangement, the condition cb^^B~^^A — c < in Lemma 
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14.21 becomes 

{Li2 — Lii + L31 — L32, L12 + I/31 — L32 — L41, 
-^13 — -^^12 + L22 — L23, L22 — L21 + L31 — L32, (26) 
L22 + L^i — L32 — L41, L31 — L41, L32 — L31 + L41 — L42} < 

where, hereinafter, {wi, . . . ,Ws} < means Wi < for i = 1, . . . , s. 

By Lemma I4.2| if ()25p and (j26p are valid, we obtain an optimal solution of the linear 
programming problem ()23p which is determined by the basic variables vetn. Hence, if 
(j26|) is valid, the corresponding q = (£1,625 £3) in and (p5|) are in Si, since in the optimal 
decomposition of q = qi + q2 + qs + ^4, qi = (0, di — 1, 1) is a vertex. 

1.2. Similarly, choosing the basic variables as veti2 = {A13, A23, A24, A31, A32, A33, A41}, 
which generates a new basic matrix B12, and from 3^2^ > 0, we obtain 

< Ai,di + d2<A2 + A3,l<A3< 2, Ai + ^2 + ^3 < 2di + ^2, 

which in turn lead to the following values for £1,62, £3 

£3 = 2, £1 = 1,..., £2 = di + ^2 - I,..., (27) 
£1 + £2 = c^i + ^2, • • • , 2(ii +d2-2. 

The condition cbi2^i2^ — c < leads to the following constraints on Lij,i = 1, . . . ,4,j = 
1,2,3, 

{L12 — Lii + L31 — L32, L12 — L32, 
-^^13 — L12 + L22 — L23, L22 — L32, (28) 
L22 — L21 + L31 — L32, L31 — L41, L32 — L42} < 0. 

1.3. Similarly, the basic variables veti3 = {A13, A23, A24, A32, A33, A41, A42} lead to 

£3 = 2, £1 = 1, . . . , (ii, 

£1 + £2 = 24 + d2 - 1, • • • , 2di + 2^2 - 2, (29) 

and 

{L12 — L42, L12 — -^^11 + L31 — L32, 
-^^13 — L12 + L22 — L23, L22 — L21 + L31 — L32, (30) 
L22 — Li2, L31 — L32 — L41 + L42, L32 — L42} < 0. 

1.4. Similarly, the basic variables veti4 = {A13, A23, A24, A33, A41, A42, A43} lead to 

£3 = 2, £1 = di + 1, . . . , 82 = di + d2 - 1, . . . , 

ei+e2 = 2di + d2,...,2di + 2d2-2. (31) 

and 

{L12 — Lii + L41 — L42, L12 — L42, 
Li3 — L12 + L22 — L23, L22 — L21 + L41 — L42, (32) 
L22 — Li2, L31 — L41, L31 — L32 — L41 + L42} < 0. 
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Let Si be the set (ei,e2,e3) defined in ([25]), (HZD, ([29]), and Tlien for r] e Si and 

an optimal sum of r] = qi + q2 + q3 + q^, since A13 = 1, qi must be the vertex (0, di — 1, 1) of 
Qi. Therefore, mm((5/i) = 2/22/1 Of course, in order for this statement to be valid, Ljj 
must satisfy constraints (j26p . (j28p . (j30p . and ()32p . We will show later that these constraints 
indeed have common solutions. 

The following three cases can be treated similarly, and we only list the conditions for 
£i,£2,£3 while the concrete requirements for are listed at the end of the proof. 



l(<5/2) = 


Ui 1 


we choose the basic variables 


vet2i = 


{Al3: 


> Ai4, 


A24, 


A3l5 


A32, A33, A41}, 


Vet22 = 


{Al3. 


: Al4, 


A24, 


A331 


A41, A42, A43}, 


Vet23 = 


{Al3: 


> Ai4, 


A24, 


A32, 


A33> A41, A43}, 


Vet24 = 


{Al3: 


> Ai4, 


A24, 


A32, 


A33> A41, A42}, 


Vet25 = 


{All: 


> A12, 


Al35 


A24, 


A31> A33, A41}, 


Vet26 = 


{Al3: 


> Ai4, 


Al5i 


A24, 


A33) A41, A43}, 


Vet27 = 


{Al2: 


. A13, 


Al4, 


Ais, 


A24, A33, A41}, 



which lead to the following elements of 52 

£3 = 1, < ei, di + d2 - I < e2,ei + 62 = di + d2 + I, ■ ■ ■ ,2di + d2 - I; 

£3 = 1, di < El, di + d2 - I < e2,ei + 62 = 2di + d2 + I, . . . ,2di + 2d2 - I; 

£3 = 1, 62 = di + d2, ■ ■ ■ , 2di + d2-l, ei + £2 = 2di + d2, ■ ■ ■ ,2di + 2^2 - 1; 

£3 = 1, £i = l,...,di, £1 + £2 = 2di + ^2, • • • , 2di + 2^2 - 1; 

£3 = 1, ei = l,...,di, 62 = d2 + l,... ,di+d2-l; 

£3 = 1, £2 = ^2 + 1, . . . , di + ^2 - 1, £2 = 2di + ^2, . . . , 2di + 2d2 - 1; 

£3 = 1, £1 = di + 1, . . . , £2 = ^2 + 1, ■ • ■ , ei + ^2 = rfi + c?2 + 2, . . . , 2(ii + (i2 - 1. 

Case 3. In order for mm(/i) = y^^, we choose the basic variables 

vet3i = {Ai5, A16, A24, A26, A33, A41, A43}, 
vet32 = {Ai3, Ai5, A23, A24, A33, A41, A43}, 
vet33 = {Ai5, A23, A24, A25, A33, A41, A43}, 
vet34 = {A12, Ai3, Ai5, A23, A24, A33, A41}, 

which lead to the following results about £1 , £2 , £3 in 53 

£3 = 1, £2 = 1,..., d2, £1 + £2 = 2di + d2, ■ ■ ■ ,2di + 2^2 - 1; 

£3 = 2, £2 = d2,... ,di +^2 - 2,£i +£2 = 2(ii + c?2 - 1,. . . ,2di + 2d2 - 2; (33) 

£3 = 2, £2 = 1,... ,(i2 - 1, £i+£2 = 2di+(i2-l,...,2di + 2d2-2; 

£3 = 2, £1 = di + 1, . . . , £2 = c?2, . . . , £1 + £2 = di + (i2 + 1, • • • , 2c?i + (i2 - 2. 
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Case 4- In order for mm(/2) = 1, we choose the following basic variables 



vet4i 


= {Aii 


Al2 


A21 


A24 


A26> A31, A41} 


Vet42 


= {Aii 


Al2 


A24 


A26 


A31) A33, A41} 


Vet43 


= {Aii 


Al2 


Al5 


A24, A26, A33, A41} 


Vet44 


= {Al2 


Al3 


A23 


A24 


A31) A33, A41} 


Vet45 


= {Al2 


A23 


A24 


A25 


A31) A33, A41} 


Vet46 


= {Al2 


A21 


A22 


A23 


A25> A31, A41} 


Vet47 


= {Al2 


A15 


A23 


A25 


A26) A33, A41} 



which correspond to the elements in S4 



£3 


= 1, 


< £1, 


< £2 


£1 + £2 = 2, . . . , ^2; 




£3 


= 1, 


£2 = 1, 


• • ,^2, 


£1 + £2 = C^2 + 1, • • • 


+ d2; 


£3 


= 1, 


£2 = 1, 


■■,d2, 


£1 + £2 = C^l + ^2 + 1, • • 


.,2di+d2- 1; 


£3 


= 2, 


£1 = 1, 


. . 


£2 = ^2, . . . , dl + (^2 - 2 




£3 


= 2, 


£2 = 1, 


. . , ^2 - 


1, £1 + £2 = C?2, . . . ,C?i 


+ ^2 - 1; 


£3 


= 2, 


£1 = 1, 


••, £2 


= 1, . . . , ^2 - 2, £1 + £2 


= 2,...,(i2-l; 


£3 


= 1, 


£2 = 1, 


. . , ^2 - 


l,£i + £2 = di + (i2 - 1 


...,2di+d2-2 



Merge all the constraints for Lij, we obtain 

Ln — L12 — L21 + -L22 < 0, 

-^13 < L23,L2i < L31 < Lii < L41, (34) 

-^22 < L12 < L32 < -^^42, -^31 = -^^32 + Lil — L42. 

The solution set for system is nonempty. For example, li = (7, —4, —5), I2 = (5, —9, 5), ^3 = 
(6, 2, 1), /4 = (8, 4, 7), which will be used for example ([I]), satisfy the conditions in 

We can also check that f = 5i U ^2 U 53 U 54 is a disjoint union for £. The lemma is 
proved. □ 

We now have the main result of this section. 

Theorem 4.4 The sparse resultant of fi, f2,Sf 1,6/2 as polynomials in y, 1/1,2/2 is not iden- 
tically zero and contains the differential resultant of fi and f2 as a factor. 

Proof. Note that ao,bo, 6ao, 6bo, which are the zero degree terms of /i, /2, Sf2 respec- 
tively, are algebraic indeterminates. As a consequence, 

Jl = {fl,f2,Sfl,Sf2) 

is a prime ideal in Q[u, y, yi,y2], where u is the set of the coefficients of /i, /2 are their first 
order derivatives. Let 

J2 = Jl nQ[u]. 
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Then J2 is also a prime ideal. We claim that 

J2 = (R) (35) 

where R is the differential resultant of /i , /2 • From c) of Theorem 12.31 R G J2 . Let T € J2 . 
Then T € Ji C [/i, /2]. From Q, the pseudo remainder of T with respect to R is zero. Also 
note that the order of T in Oj, 6j is less than or equal to 1. From a) and b) of Theorem 12.31 
R must be a factor of T, which proves (|35p . 

From Lemma 14.31 the main monomials for fi, f2,Sfi,^f2 are the same as those used to 
construct 5i , ^2 , ^3 , ^4 in ([5]) . As a consequence, we have Si C Si. For q & Si\ Si, q must 
be in some Si, say q G 5*2. Then from Lemma 14.31 the monomials in (M(g)/mm((5/2))5/2 
are contained in £. By Corollary 13.31 the sparse resultant matrix of /i, /2, Sfi, 6/2 obtained 
after move q from ^2 to Si is still nonsingular. Doing such movements repeatedly will lead to 
Si = Si,S2 = S2,S3 = S3, Si = S4. As a consequence, the sparse resultant is not identically 
zero. 

From (jlSp . we have TZ £ Ji which implies 7^ € J2. Since R is irreducible, R must be a 
factor of TZ. □ 

4.3 Example ([1]) revisited 

We show how to construct a nonsingular algebraic sparse resultant matrix of the system 
{gi,g2,6gi,6g2}, where gi,g2 are from 

Using the algorithm for sparse resultant in [21 [3], we choose perturbed vector 6 = 
(0.01,0.01,0.01) and the lifting functions h = (7, -4,-5), ^2 = (5,-9,5),Z3 = (6,2, l),/4 = 
(8,4,7), where k corresponds to Qi defined in ([22]) with di = 6,2 = 2. These lift functions 
satisfy the conditions ([Mj) . 

By Lemma 14.31 the main monomials for gi,g2,dgi,Sg2 are identical with those given in 
Section [3.31 Let Si, S2, S3, S4 be those constructed as in the proof of Lemma l4.3i After the 
following changes 

move {y2yiy^,y2yiy'^} in S3 to 5i, 
move {y2y^,yiy^,y^,y^} in 5*4 to S3, 
move {y2yiy,y2yi} in 54 to Si, 

we have Si = Si, i = 1, ... ,4. Then by Corollary 13.31 the sparse resultant matrix constructed 
with the original 5i , 5*2 , 5*3 , 5*4 is nonsingular and contains the differential resultant as a 
factor. 

5 Conclusion and discussion 

In this paper, a matrix representation for two first order nonlinear generic ordinary differen- 
tial polynomials /i,/2 is given. That is, a non-singular matrix is constructed such that its 
determinant contains the differential resultant as a factor. The constructed matrix is further 
shown to be an algebraic sparse matrix of /i, /2, 5/2 when certain special lift functions 
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are used. Combining the two results, we show that the sparse resultant of /i, /2, Sfi, 8/2 is 
not zero and conatins the differential resultant of /i, /2 as a factor. 

It can be seen that to give a matrix representation for n + 1 generic polynomials in n 
variables is far from solved, even in the case of n = 1. Based on what is proved in this paper, 
we propose the following conjecture. 

Conjecture. Let V = {/i, /i, . . . , fn+i} be n + 1 generic differential polynomials in n 
indeterminates, ord(/j) = Sj, and s = X^"=o^«- Then the sparse resultant of the algebraic 
polynomial system 

h,6h,...6'-''^h,...,fn+l,6fn+l,... 6'-'" fn+1 (36) 

is not zero and contains the differential resultant of "P as a factor. 
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